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Abstract 



The Hamiltonian model of quantum error correction code in the literature is often 
constructed with the help of its stabilizer formalism. But there have been many 
known examples of nonadditive codes which are beyond the standard quantum error 
correction theory using the stabilizer formalism. In this paper, we suggest the other 
type of Hamiltonian formalism for quantum error correction code without involving 
the stabilizer formalism, and explain it by studying the Shor nine-qubit code and 
its generalization. In this Hamiltonian formulation, the unitary evolution operator 
at a specific time is a unitary basis transformation matrix from the product basis 
to the quantum error correction code. This basis transformation matrix acts as an 
entangling quantum operator transforming a separate state to an entangled one, 
and hence the entanglement nature of the quantum error correction code can be 
explicitly shown up. Furthermore, as it forms a unitary representation of the Artin 
braid group, the quantum error correction code can be described by a braiding 
operator. Moreover, as the unitary evolution operator is a solution of the quantum 
Yang-Baxter equation, the corresponding Hamiltonian model can be explained as 
an integrable model in the Yang-Baxter theory. On the other hand, we generalize 
the Shor nine-qubit code and articulate a topic called quantum error correction 
codes using Greenberger-Horne-Zeilinger states to yield new nonadditive codes and 
channel-adapted codes. 
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1 Introduction 



Quantum information and computation [JJ [2] is a new interdisciplinary field combining 
(quantum) physics, (advanced) mathematics and (modern) computer science, and one 
of its great aims is to build a real quantum computer on which quantum algorithms 
can successfully run. Quantum error correction codes are exploited to protect quantum 
information from various kinds of noise and perform a large-scale quantum computa- 
tion with imperfect quantum gates. They were originally invented by Shor, Steane, 
Calderbank, et al. [31 01 El [6] , for examples, the Shor nine-qubit code [3], the Steane 
seven-qubit code and the perfect five-qubit code saturating the quantum Ham- 

ming bound 0IB]- Additive quantum codes analogous to classical linear codes can be 
described in the stabilizer formalism which have two well known forms [HlEl] and [11] . 
Quantum error correction conditions [7J [12] play a fundamental role in the quantum 
error correction theory as a good guidance of devising good quantum error correction 
codes. Fault-tolerant quantum computations [El [HI [El [161 [T7] study how to per- 
form reliable quantum computation with the help of quantum error correction codes, 
and the threshold theorem [TBI E3 can De used to evaluate various proposals for fault- 
tolerant quantum computations and guide physicists to devise reasonable experiments 
in quantum information and computation. 

Quantum computer is a physical system described by quantum mechanics, i.e., its 
dynamics determined by the Hamiltonian in the Shrodinger equation. The focus of the 
present paper is to study the Hamiltonian model underlying quantum error correction 
code. In Kitaev's toric code [14] , the Hamiltonian is a linear combination of elements 
of the stabilizer group, its degenerate ground state represents a quantum error correc- 
tion code, and there is a gap between its ground state and its first excited state to 
protect this code from environment noise. Subsequent to the toric code is topological 
quantum computing |14[ [20] in which qubits are anyons (quasiparticles obeying the 
braid statistics) and quantum gates form unitary braid representations. In condensed 
matter physics, Abelian or non-Abelian anyons can be created in the fractional quan- 
tum Hall effect [21] • Recently, encoding quantum information into a subsystem of a 
physical system has been found to be a most general method for protecting quantum 
information from decoherence, and it is summarized in operator quantum error correc- 
tion [12] |22~1 123] . Similar to the Hamiltonian model of the toric code, two examples 
for how to construct a Hamiltonian model for operator quantum error correction sub- 
system codes are presented by Bacon [24] . Remarkably, the Hamiltonian model of the 
toric code [25] on a four dimensional lattice and the Hamiltonian model of the sub- 
system code [24] on the three-dimensional cubic lattice may be theoretical candidates 
for self-correcting quantum memory where the robust storage of quantum information 
is guaranteed by physical properties of the system (i.e., the Hamiltonian and bound- 
ary conditions or topology). Besides, Hamiltonian models for self-correcting quantum 
memory, fault-tolerant quantum computation can be described by a physical system 
with a time-dependent Hamiltonian [26 ] 127 ] [28] . In this paper, however, we will discuss 
quantum error correction code in a different Hamiltonian formalism from what have 
been introduced as above. We study quantum states themselves which represent a 
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quantum error correction code, rewrite them by a unitary transformation on the prod- 
uct basis, and then derive the Hamiltonian to determine this unitary transformation. 
The stabilizer formalism is not involved in our time-independent Hamiltonian formu- 
lation. We will explain our motivation in detail with simple examples in the following 
sections. 

In the literature, there have been many new kinds of quantum error correction 
approaches, for examples: nonadditive codes [23 EH1 ED E21 E3] beyond the stabilizer 
formalism, entanglement-assisted quantum error correction [M], topological color codes 
[33] , subsystem codes [M1EZ]) channel-adapted codes [3H1 E3 HO] , etc. Creating new 
quantum error correction code is not a primary goal of this paper, whereas revisiting 
quantum error correction codes from a new point of view is our first point. Quantum 
error correction code is represented by a unitary basis transformation matrix from the 
product basis to the entangling basis in terms of quantum error correction code, and 
this explicitly shows the entanglement nature of quantum error correction codes (i.e., 
"fight entanglement with entanglement" by Preskill [15\ 116]). A possible experiment 
realization of such a basis transformation matrix can be referred to [H] . We are only 
concerned about quantum error correction code itself instead of its underlying physical, 
informational, algebraic, graphical, topological, and geometrical mechanism. On the 
other hand, we will indeed articulate a new topic called "quantum error correction 
codes using Greenberger-Horne-Zeilinger (GHZ) states" which include as examples the 
Shor nine-qubit code [3], amplitude damping channel-adapted codes [381 E3 HQ] , and 
nonadditive codes \29\ HQ] . It is possible to invent a theoretical framework for devising 
quantum error correction codes using GHZ states based on recognizable properties of 
GHZ states. GHZ states [42 1 1431 HI] are usually assumed to be maximally entangled 
states in various entanglement measure theories [TJ [2] . 

As a remark, our research represents a further development of a recent study pre- 
sented infMlHElHSHHlHnilSOliniEaEaiMlEa We exploit observations and 
techniques in these articles (especially [52] ) to study quantum error correction, which 
has not been done before to the best of our knowledge. They mainly study interdis- 
ciplinary themes arising from quantum information and computation, low dimensional 
topology [57], the Yang-Baxter equation EHl ED] and (almost-complex) differential 
geometry. For examples: connections among quantum entanglements, topological en- 
tanglements and geometric entanglements as well as integrable quantum computation 
are discussed in [451 "l6l Wl\ 11*8] ; a possible link [52] between quantum error correction 
and topological quantum computing can be set up via the Bell matrix, which forms a 
unitary braid representation and acts as a unitary basis transformation matrix from 
the product basis to GHZ states. In the sense of quantum information and computa- 
tion, these papers articulate integrable quantum computation as well as topological-like 
(partial topological) quantum computation [5H ESI EE] : the former is a new approach 
for quantum computation using integrable models [27J H3 EH1 ESI ED] constructed via 
solutions of the quantum Yang-Baxter equation (the Yang-Baxter equation with the 
spectral parameter); and the latter exploits low dimensional topology [I5T] , unitary braid 
gates and non-braid gates. In our Hamiltonian formulation, Shor's code (and its gener- 
alization) can be explained as a braiding operator, and the corresponding Hamiltonian 
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model can be recognized as an integrable model. 

We hereby summarize our main result which is new to our knowledge. 1) Quantum 
error correction code is described in a new Hamiltonian formulation, and the unitary 
evolution operator at the specific time is a unitary basis transformation from the prod- 
uct basis to the quantum error correction code. 2) In our Hamiltonian formulation, 
Shor's code and its generalization can be recast as braiding operators. 3) Quantum er- 
ror correction codes using GHZ states are analyzed as an independent topic for yielding 
new nonadditive codes or channel-adapted codes. The plan of this paper is organized as 
follows. In Section 2, with the helpful formalism of GHZ states, Shor's nine-qubit code 
is described in the Hamiltonian formulation and then explained in the language of the 
braid group. In Section 3, we discuss quantum error correction codes using GHZ states, 
study generalized Shor's codes as examples, and make comments on nonadditive codes 
and channel-adapted codes. Finally, we conclude this paper by presenting interesting 
problems for further research. 

2 Notations, motivations and Shor's nine-qubit code 

Notations for most symbols in this paper are introduced, and an overview is made 
on physical, informational and mathematical properties of GHZ states. Hamiltonian 
formulations of Shor's nine-qubit code are explored in detail, so that our motivations 
for the Hamiltonian formulation of quantum error correction codes (without exploiting 
the stabilizer formalism) can be understood well. Besides, the original motivation takes 
root in our previous work [52| recasting GHZ states as unitary solutions of the Yang- 
Baxter equation (or the braid group relation). 

2.1 Notations 

In quantum information and computation [TJ [2] , a two dimensional Hilbert space H.2 
over the complex field C is called a qubit, for example, a|0) + /3\1), a, (3 £ C where |0) 
and |1) form an orthonormal basis of TL2 — C 2 and are usually chosen as eigenvectors 
of the Pauli operator Z: Z\0) = |0) and Z\l) = — 11). The Pauli matrices have the 
conventional form, 

*-(;;)• *-(i-0- r = zx ={-i 1) « 

which denote the bit-flip, phase-flip, and bit-phase flip operation on a qubit. 

The symbol 11 2 denotes a 2-dimensional identity operator or 2 <g> 2 identity matrix. 
The notations A® n and |a)® n denote the following ra-fold tensor products, 

A® n = A®---A, |a)® n = \a) <g) ■■■ <g> |a) . (2) 

" * ' V v ' 

n n 

An n-fold tensor product, in terms of the identity operator I2 and Pauli matrices 
X, Y, Z, has a simpler notation in which the lower index of the Pauli matrix labels its 
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position in this tensor product, for example, Z\Z 2 and X^X^X^X-jX^X^ in a 9-fold 
tensor product describing the form, respectively, 

Z 1 Z 2 = Z®Z®{\ 2 )® 7 , X 4 X 5 X 6 X 7 X 8 X 9 = (l 2 f 3 ®X® 6 . (3) 

The symbols M and B are exploited in the entire paper. The M is an anti-Hermitian 
operator = —M satisfying M 2 = —Id with Id denoting the identity operator, and 
B = ei M . But M and B will have different presentations depending on where they 
appear. If M is a /s-fold tensor product in terms of Pauli matrices X, Y, Z, then we 
introduce a notation for an n-fold tensor product M, by 

Mi = (l2) 0i_1 ®M® (l 2 ) n ~ k ~ i+ \ l<i<n+l-k (4) 

where the lower index of can be relabeled for convenience. In [51} 152] , the symbol 
M is called the almost- complex structure in the (almost-complex) differential geome- 
try, and it can yield an anti-Hermitian representation of the extraspecial two-groups. 
As B = e* M is a solution of the Yang-Baxter equation [581 [59] and can generate a 
unitary braid representation 7r n of the Artin braid group B n [57], it is called the Bell 
matrix, see [ITJ EE] EH EU [52]. Since unitary braids are used as quantum gates in 
topological quantum computing [HI [20] and extraspecial two-groups are exploited in 
quantum error correction [10] , the formula B = e^ M suggests there may exist a link 
between quantum error correction and topological quantum computing, which is the 
main proposal of our previous paper [52] . 

2.2 An overview on Greenberger-Horne-Zeilinger states 

The GHZ states [42 [ I43 [ [4"4"] are a multipartite generalization of the bipartite maximally 
entangled Bell states, and are defined by various known entanglement measures as max- 
imally entangled. They play crucial roles in both fundamental problems and practical 
applications of quantum information theory. The Bell theorem [f"T] for incompatibility 
between quantum theory and classical deterministic local models is expressed in the 
form of inequalities (Bell inequalities) among various statistical correlations, whereas 
the GHZ theorem [421 [4*3] asserts that the quantum correlations represented by the 
GHZ states allows us to describe the Bell theorem in terms of equalities and to test the 
expected incompatibility only by perfect correlations. The GHZ states are the simplest 
multipartite maximally entanglement sources, and have been widely exploited in the 
study of quantum information and computation, for examples, multiparty quantum key 
distributions [62] and quantum teleportation [63j . Besides, GHZ states are often called 
cat states acting as ancillas in fault-tolerant quantum computation |13[ [T5l [T6l [17] . 

Besides the fundamental importance of GHZ states in quantum physics and infor- 
mation, they have been found to posses (beautiful and deep) topological, algebraic and 
geometric properties. Aravind [45] observed that there are similarities between the GHZ 
states and knot configurations [57], by identifying the measurement of a specific state 
of a particle with cutting the corresponding link component. Furthermore, we describe 
GHZ states by the Bell matrix B = e^ M which has become a common topic among 
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quantum error correction, topological quantum computing, the braid group, the Yang- 
Baxter equation, and (almost-complex) differential geometry jl6j WI \ HHj H9 | 1501 |5T | [52] . 
The GHZ states of n-qubit are defined to have the form, 

-^=(|s) ± s = ■ ■ ■ i n , s = i 1 i 2 ---i n (5) 

where ij = 0, 1, j = 1, • • • , n and ij + i j = 1 with the Abelian addition modulo 2. In 
[52] , the GHZ states of the n-qubit are generated by the exponential function B = e^ M 
acting on the product basis, for example, 

-^(|0)® n + |l)® n ) = B\0)® n , -^L(|O) 0n - |l)® n ) = J B" 1 |0)® n , (6) 
V2 v 2 

where M = -Y®X® n and = e~^ M denotes the inverse of the B matrix. Note that 
B = e^ M is called the Bell matrix only if it is a unitary solution of the Yang-Baxter 
equation |58[ [5TJ] (or forms a unitary braid representation .57] ) . 

2.3 Hamiltonian formulation of Shor's nine-qubit code 

The first quantum error correction code was found by Shor [3] , and it describes a logical 
qubit state in terms of nine-qubit states in the way 

|0) -» \0)l = -A=(|000) + |111>) ® (|000) + |111» &> (| 000) + |111>), 
2y 2 

|1) - = i(|000) - |111» (|000) - |111» ® (|000) - |111>) (7) 

where two GHZ states -^(|000) ± 1 1 1 1 ) ) are repeatedly used and hence Shor's code is 
also called the repetition code. In the stabilizer formalism for additive quantum 
error correction codes, the encoded logical qubit a\0) l a, (3 £ C is uniquely de- 

termined as a common eigenvector with the eigenvalue 1 of the following eight stabilizer 
operators, 

Z1Z2, Z2Z3, Z4Z5, Z$Zq, ZfZg, Z$Zg, X1X2X3X4X5XQ, X^X^XqXjXsXq (8) 

which form an Abelian group called the stabilizer group. The logical bit-flip operation 
X and logical phase-flip operation Z usually have the form 

X = Z® 9 , ~Z = X m . (9) 

In the following, we firstly describe Shor's nine-qubit code in our Hamiltonian for- 
mulation and then present an interpretation in the language of the braid group. Denote 
the symbol M by M = —Y®Y®Y in this subsection. It is a transition operator between 
|000) and | 111) , 

M|000) = I 111) , Af|lll) = -|000) (10) 
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due to — Y\0) = |1) and — = — 10), and it is an anti-Hermitian operator M' = —M 
determining e* M to be a unitary operator, 

e? M = -^(if 3 + M), M 2 = -lf. (11) 
V2 

Furthermore, introduce Mi, M4, M7 as 

Mi = M®lf 6 , M 4 = lf 3 ®M<g>lf 3 , M 7 = if 6 <g) M (12) 

and denote the summation of Mi, M4, M7 by M t = Mi + M4 + M7. After some simple 
algebra, the Shor nine-qubit code has a new compact formulation, 

|0> L = e ? M '|0>® 9 , = e -f A/ *|0)® 9 (13) 

which leads to the bit-phase flip operation Y on the encoded qubit, 

|0) L = ef M *|l) L , |l) L = e -f M *|0) L . (14) 

Namely, Y has the form 

Y = e% Mt = MxM^M-j = —Y® 9 (15) 

where the minus sign in the front of Y® 9 partly explains the reason why in the literature 
the logical bit-flip operator X and the logical phase-flip operator Z have the formulation 
© because of Y = ZX and XZ = -ZX. 

Now let us associate the formulation (I14p of the Shor nine-qubit code with the 
Hamiltonian defined by H = iM t . A logical qubit basis, |0)l and is explained as a 
unitary evolution of the product basis state |0)® 9 (or |1)® 9 ), with the unitary evolution 
operator U(8) = e~ ldH determined by the Hamiltonian H, namely, 

|0>L = t/(^)|0r, \1) L = U(-^)\0)® 9 = -U(l)\l)*> (16) 
which leads that the encoded qubit a|0)x + (3\1)l has a form 

a\0) L +[3\l)L = U^)(a\0)® 9 -f3\l)® 9 ), a,/3eC. (17) 
The Shrodinger equation has the form with the unitary evolutional solution U(6), 

i^{0) = H^{6), H = -i{Y x Y 2 Y3 + Y 4 Y 5 Y 6 + Y 7 Y 8 Y 9 ) (18) 

where the Planck constant h = 1, 9 is regarded as the time variable and ip(9) represents 
the wave function in quantum mechanics. The three-body spatially local Hamiltonian H 
can be simulated on various lattice models, but all of which may be equivalent to a one- 
dimensional spin chain because only the Y operation is involved in this Hamiltonian. 
Note that we do not impose any boundary conditions on the Shrodinger equation (|18j) . 
mainly because they are determined by which type of lattice model to be chosen. 
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In our Hamiltonian formulation, obviously, Shor's nine-qubit code is not an eigen- 
vector of the Hamiltonian H and is not its ground state. This is essentially different 
from Hamiltonian models [14] [24] respectively for the toric code and subsystem code, 
which are constructed with the help of the stabilizer formalism. In spite of this fact, we 
expect our Hamiltonian formulation to be helpful for exploring the topics of how to en- 
code and decode Shor's code, how to perform the fault-tolerant quantum computation 
with Shor's code, and especially how to construct Hamiltonian model of nonadditive 
code without using the standard stabilizer formalism. In the next subsection, further- 
more, we will recast the Shor nine-qubit code as a braiding operator, and explain the 
motivation for our Hamiltonian formulation from the different point of view. 

2.4 Artin braid group, Yang— Baxter equation and Shor's code 

The Artin braid group B n [57] has the generators cij, i = 1, ■ ■ ■ ,n — 1 satisfying the 
braid group relations 1) and commutative relations 2): 

1) OiOi+\Oi = a i+ i(Jicri, 2) ^aj = ajdi, \i - j\ > 1. (19) 

In terms of the identity matrix 11 2 and a /c-fold tensor product M involving Pauli matri- 
ces, a unitary braid representation n n of the Artin braid group B n can be constructed 
in the way, 

B i = ir n (a i ) = (\ 2 )® i - 1 ®e^ M ®{\ 2 )® n - k - i+ \ \<i< n + l-k. (20) 

where e< M satisfies the Yang-Baxter equation [581 [59| [52] . 

(G Id)(Id ®G)(G® Id) = (Id ® G)(G ® Id)(Id <g> G), (21) 

where G = e* M and Id = I2 and which is a presentation of the braid group relation 
<jj<7j_|_icrj = (Tj+icrjcrj, see [S2] for more details. Besides, the quantum Yang-Baxter 
equation is defined as the Yang-Baxter equation with the spectral parameter, 

(G(x) ® Id)(Id <g> G{xy)){G(y) <g) Id) = (Id ® G(y))(G(xy) ® Id)(Id ® G(x)), (22) 

where x, y are the spectral parameter and which can be referred to our previous papers 
|47[ [4"5"1 [5"0l [5"T| [5"2"] . It is well known that in the literature an integrable model can be 
constructed using a solution of the quantum Yang-Baxter equation, see [60] . 

In Figure 1, the diagram without boxes around crossings is a diagrammatic rep- 
resentation of the Yang-Baxter equation pip or (|22H in which different strands are 
allowed to represent Hilbert spaces of different dimensions, for example, thin strands 
acting on H2 (i.e., a qubit) and thick strands acting on the tensor product of TL2 (i.e., at 
least two qubits). On the other hand, G = e^ M can be explained as a universal quan- 
tum gate |64[ [4*6] which transforms a separate state to a maximally entangled state, for 
example, a GHZ state. The diagram with boxes is hence regarded as a diagrammatical 
identity between two quantum circuits where universal quantum gates are labeled by 
boxes with inputs and outputs. In other words, we are presenting two sorts of inter- 
pretations on the same diagrammatical object: the one is the Yang-Baxter equation 
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[G®Id) (Id ® G) [Id <g> G) = (/d ® G) (G <g) 7d) (/d <g> G) 



Figure 1: Yang-Baxter equation and quantum circuits 



(or the braid group relation) [571 ESI E2 E0] and the other is quantum information 
and computation [H[2]. We are trying to build a bridge between quantum information 
theory and Yang-Baxter equation (or low dimensional topology), which is the under- 
lying motivation of a series of papers [UJ HHJ EQl ED E21 E31 EH ESI [56] . For examples, 
we have recognized Werner states [HS] as a sort of rational solutions of the quantum 
Yang-Baxter equation [53]; in the formulation of GHZ states [511 152] . found mathe- 
matical structures including unitary braid representations, extraspecial two-groups and 
almost-complex structure; described quantum teleportantion protocol [63j by a braiding 
operator [541 155]; and recast diagrammatical quantum information protocols involving 
maximally bipartite entangled states as the extended Temperley-Lieb diagrammatical 
configuration [541 [551 [56] . In the following, we refine the braid group relation from our 
Hamiltonian formulation of the Shor nine-qubit code. 

Choose M to be the form M = —Y (g) X <g> X and construct the unitary braid 
representations B\,B±,B 7 by 

B x = e\ M ® If 6 , B 4 = If 3 ® e* M (8) if 3 , B 7 = if 6 ® e 5 M , (23) 

and then rewrite Shor's nine-qubit code into a compact formulation 

\0) L = B^B-rlO)® 9 , \1) L = B^B^B^IO)® 9 . (24) 

The B\,B±,Bi are presentations of the braid group generators, and have a diagram- 
matic representation, see Figure 2. Every straight strand denotes a two-dimensional 
Hilbert space H2 (i.e., a qubit), and the crossing denotes e* and acts on Tif 3 . For 
convenience, we can assign TL2 to the first strand of the crossing and TL® 2 to the second 
strand. Therefore, we read the Shor nine-qubit code from the braiding configuration 
as a product of three braids Bi,B±,Bj. 

On the other hand, in terms of Mt = M\ + M4 + Mj where Mi , M4, Mj have a form 
using the convention ([S]) , 

M 1 = -Y 1 X 2 X 3 , M 4 = -Y 4 X 5 X 6 , M 7 = -Y 7 X 8 X 9 , (25) 

we obtain the other Hamiltonian formulation of Shor's code similar to (|14j) but the 
logical bit-phase flip operation Y given by 

Y = -YiX 2 X 3 Y4X 5 X 6 Y 7 X 8 X 9 . (26) 
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Figure 2: Artin braid group generators for Shor's nine-qubit code 



In view of our previous work on integrable quantum computation [47} [4"51 [5"U1 [5"T| [52"] , the 
unitary evolution operator U(9) = e~ l8H with H = iM% is a solution of the quantum 
Yang-Baxter equation (|22p with the spectral parameter 9, and hence the corresponding 
Hamiltonian model is an integrable model determined by the solution of the quantum 
Yang-Baxter equation [58J EH EQ] . 

As a remark about our unitary braiding description of Shor's nine-qubit code, it 
provides an interesting example for the observation [52] that there exists a possible link 
between quantum error correction and topological quantum computing via the unitary 
braid representation using anti-Hermitian representation of extraspecial two-groups. 
Besides, it is possible to make connections between our work and a geometric paradigm 
of Shor's code [66] . 

3 Quantum error correction codes using GHZ states 

Quantum error correction codes using GHZ states are defined as those codes in terms 
of GHZ states or tensor products of GHZ states, while generalized Shor's codes using 
GHZ states are specified to be repetition codes as a straightforward generalization of the 
Shor nine- qubit code, i.e., only involving GHZ states of the type, l/\/2(|0)® n ± \l)® n ). 
They are degenerate quantum error correction codes in which different errors may be 
corrected by the same correction operation. Besides the Shor nine-qubit code, some 
nonadditive codes [201 HO] and channel-adapted codes [351 E3 SO] are good examples 
for quantum error correction codes using GHZ states. The motivation for our articula- 
tion of this topic is based on the observation that GHZ states have nontrivial physical, 
informational, and mathematical properties (see our overview on GHZ states in Sub- 
section 2.2), and these properties are expected to be very helpful in the construction 
of interesting quantum error correction codes using GHZ states. 

In the literature, the triple [n, k, d] of natural numbers represents a quantum error 
correction code [I]. It is a unitary mapping Q from Ti.f k to H.® 11 (k < re), and encodes 
information of &-qubit quantum states into n-qubit quantum states in order to detect 
and correct t-qubit errors, where t = for d even and t = for d odd. Quantum 
states in this 2 fc -dimensional subspace Q7if k are called codewords. For example, Shor's 
nine-qubit code is usually denoted by [9, 1, 3] and it can correct any single-qubit errors, 
while the logical qubit a|0)i + ot,(3 £ C represents a codeword in Shor's code. 

In this section, for simplicity, we denote quantum error correction codes [n, k, d] by 
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[n, k], and focus on whether or not they are able to correct arbitrary single-qubit errors, 
i.e., t = 1. For example, the Shor code is denoted by [9,1]. Note that we emphasize 
our perspectives of constructing quantum error correction codes using GHZ states and 
describe them in our Hamiltonian formulation, but leave the systematic construction of 
procedures for encoding, error syndrome measurement, error correction and decoding 
of these codes elsewhere. 

3.1 Generalized Shor's codes using GHZ states 

First of all, we introduce symbols and names necessary for constructing generalized 
Shor's codes using GHZ states. The type of a GHZ state ± \s)) © is specified by 

the number of its qubits and the sign ±. Generalized Shor's codes exploit different types 
of GHZ states to encode one qubit into codewords of iV-qubit. The same GHZ states 
form a block in the codeword, and different blocks correspond to orthogonal subspaces 
of Hf N . The symbol / is the number counting types of GHZ states (or the number 
of blocks) in the codeword; i denotes the i-th type of GHZ state (or i-th block) in the 
codeword, i = 1, q% counts the number of qubits in the i-th type of GHZ state; n{ 

is the repetition number of the i-th type of GHZ state in the i-th block of the codeword; 
B = Y^,i=i n i counts the number of GHZ states in a codeword; N = Yli=i n iQi ^ s the 
total number of qubits in the codeword. 

For example, the Shor nine-qubit code exploits GHZ states of three-qubit and re- 
peatedly use them three times in a codeword: I = 1, q\ = 3, n\ = 3, B = 3 and N = 9. 
On the other hand, B = 3 is the number counting all phase-flip single-qubit errors 
(i.e., Z-errors), and 2N = 18 counts all bit-flip and phase-flip errors (i.e., X-errors and 
Y-errors). Hence Shor's code satisfies the quantum Hamming bound condition, 

2*(1 + B + 2N) < 2 N (27) 

where t = 1, and which suggests the Shor code as a degenerate code and able to correct 
some two-qubit errors. Now let us discuss an application of the above bound inequality 
at t = 1. Obviously, the lowest bound for N has to be N > 5. 

At N = 5, we can either construct a logical qubit using GHZ states of five-qubit, 

|0) L = -J=(|0)« 5 + |l)n \1) L = ±(\0)®* - (28) 

where B = 1, n\ = 1 and q\ = 5, or have a logical qubit in terms of Bell states (i.e., 
GHZ states of two-qubit) and GHZ states of three-qubit, 

|0) L = i(|0)® 2 + |1) 02 )(|O) 03 + |1) 03 ), |1) L = ^(|O) 02 - \l) m ){\Q) m - \l) m ) (29) 

where B = 2, m = ri2 = 1, qi = 1 and q2 = 1. 

At N = 6, we can construct a logical qubit in terms of GHZ states of six-qubit (or 
four-qubit or three-qubit) or Bell states, for examples, 

\0) L = ^(|0)® 2 + |1)® 2 )(|O) 02 + |1)® 2 )(|0)® 2 + |1)® 2 ), 
2v 2 
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\1) L = -^(|O> 02 - |1> 02 )(|O> 02 - |1> 02 )(|O> 02 - |1) 02 ). (30) 

where B = 3, n\ = 3 and f/i = 3. 

In view of our Hamiltonian formulation of the Shor nine-qubit code, generalized 
Shor's codes have the following Hamiltonian formulation, 

|0) L = et^ M "|0)^, \l) L = e-5&=i M P\0)* N (31) 

where the Hamiltonian = y/—lM^ determines how to yield the i-th type of GHZ 
state in the i-th block of the codeword and has the form 

rii 

#(0 = V^J^Mj . (32) 
i=i 

Different M- are commutative with each other because the corresponding subspace 

of each type of GHZ state in the codeword is orthogonal with those for other types. 

(i) 

Note that the lower indices of M- are different from that exploited in the Hamiltonian 
formations of the Shor nine-qubit code in Subsections 2.3 and 2.4. 

Note on a link between the cat code [re, 1] and quantum error correction codes using 
GHZ states. The cat code [n, 1] has the form |0} L = \0)® n and \1) L = |l) 0n leading to 

-L(|0)z ± |1) L ) = e ±i M \0) L , |0)x + = M(\0)l ~ (33) 

in which M can take a form M = —Y <£> X® n or other possible forms. Besides, it 
is interesting to compare generalized Shor's codes using GHZ states with generalized 
Shor's subsystem codes [36], because they are constructed under completely different 
motivations and terminologies. 



3.2 Comments on channel-adapted codes and nonadditive codes 

As has been shown in Subsection 2.2, GHZ states are maximally entangled multipartite 
quantum states and posses (simple and deep) algebraic, topological and geometric 
properties. These properties have been used to explore connections among quantum 
information and computation, the Yang-Baxter equation, low dimensional topology 
and differential geometry, see [52] . On the other hand, GHZ states have been exploited 
to yield quantum error correction codes [2SJ EH1 [Ml HO] • Some of them [3E1 ES] can be 
described in the stabilizer formalism, i.e., additive codes, but many of them [29t 140] are 
nonadditive codes which do not have classical analogues and do not satisfy the quantum 
Hamming bound (|27|) . Note that nonadditive codes using GHZ states are called self- 
complementary codes in [291 HO]- Some codes using GHZ states (either additive or 
nonadditive) are channel-adapted quantum error correction for the amplitude damping 
channel [38 ] 1391 [4*0] . However, it remains an open problem how to create nonadditive 
codes using GHZ states in a unified theoretical framework, though interesting examples 
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have been found. We expect that algebraic (or topological or geometrical) properties 
of GHZ states [52] are helpful for solving this problem. This is the motivation for our 
articulating quantum error correction codes using GHZ states. 

As examples, we analyze two codes [38j ES] devised for the amplitude damping 
channel in our Hamiltonian formalism. They may be helpful for seeking a theoretical 
framework underlying nonadditive codes using GHZ states. The first one [38] is the 
code [4, 1] having a logical qubit spanned by 

|0)l = ^=(10000)) + |1111)), \1) L = ^=(10011) + |1100)). (34) 

This code can be rewritten as the unitary basis transformation e _i * H , 

) L ,\1)l) =e-^ H (|0000),|0011)) (35) 



where iH = Y ® X® 3 is called the almost-complex structure and the Bell matrix 
B = e~' l 4 H forms a unitary braid representation, see [52]. Note that the Bell matrix B 
can generate all the GHZ states of four-qubit from the product basis, i.e., it has all the 
information of GHZ states of four-qubit. In this sense, the Hamiltonian formulation 
of the code [4, 1] suggests: it appears to work only with two GHZ states of four-qubit, 
but in fact it involves all other GHZ states of four-qubit through the Bell matrix B. 
This may be one of the reasons why it violates the quantum Hamming bound (I27p but 
works well for the amplitude-damping channel. The second example is the code [6, 2] 
which has the form in the Hamiltonian formulation, 

(|00) L , \01) L , |10)i, |ll)i) = e - ^(|000000), (001001), (000110), (110000)), (36) 

where H = —iY<giX® 5 . A similar analysis can be made: the unitary evolution operator 
e~ l i H has the information of all GHZ states of six-qubit, and hence the code [6, 2] 
actually encodes two-qubit information into the entire Hilbert space spanned by 2 6 = 
64 GHZ basis states instead of its four-dimensional subspace. In this sense, we may 
understand why a 12-dimensional subspace can be encoded into the GHZ states of 
8-qubit and protected by a nonadditive code [4*0] . 

Furthermore, we discuss how to protect a one-dimensional space by quantum error 
correction codes using GHZ states. The quantum Hamming bound condition (|27p at 
t = has a perfect solution in terms of B = 1, N = 3 saturating the bound, i.e., 
(l + l + 2x3)<2 3 . Let us encode the one-dimensional space using a GHZ state 
of three-qubit, 

\t/>) L = 4= (I 000 ) + I 111 )) = e~^ H \0) m , H = -iY®Y®Y. (37) 
v2 

There are three bit-flip errors, one phase-flip error and three bit-phase flip errors giving 
rise to the other seven GHZ states of three-qubit, 

X 1 \^) L ,X 2 \^) L ,X 3 \^) L ,Z 1 \iP) L ,Z 1 X 1 \i;) L ,Z 2 X 2 \iP} L ,Z 3 X 3 \tl;) L , (38) 
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which form an orthonormal basis of Ttf together with They can be recast in our 

Hamiltonian formulation 

(Id, X 1 ,X 2 ,X 3 ,Z 1 , Z 1 X 1 , Z 2 X 2 , Z 3 X 3 )\ip} L 

= e"^ H (|000), |011), 1 110) , -|111), -|100), -|010), — 1001>) (39) 

where four minus signs can be absorbed by rescaling the phase-flip Z-error. Besides, 
we can exploit the standard procedures of encoding, error syndrome measurement, 
correcting and decoding devised for the stabilizer codes [TTj . 

4 Concluding remarks and outlook 

Motivated by the problem how to construct a Hamiltonian model for nonadditive quan- 
tum error correction code, we suggest a new Hamiltonian formulation of quantum error 
correction code without appealing to its stabilizer formalism (if it exists). As a remark, 
in Kitaev's fault tolerant schemes |14| . thermal effect may not be easily surmountable 
H which can destroy the expected fault-tolerance, see [67], and hence seeking for new 
Hamiltonian models of quantum error correction codes remains a fundamental problem 
in the current research of quantum information and computation. 

In our Hamiltonian formulation, the unitary evolution operator at a specific time 
is a unitary basis transformation matrix from the product basis to the quantum error 
correction code, and we explain this by studying examples including Shor's nine-qubit 
code and its generalization. Remarkably, as this basis transformation matrix is a solu- 
tion of the Yang-Baxter equation, the quantum error correction code can be explained 
as a braiding operator, and the Hamiltonian model is an integrable model determined 
by the solution of the quantum Yang-Baxter equation. On the other hand, we ar- 
ticulate the topic called quantum error correction codes using GHZ states, in which 
new nonadditive codes and channel-adapted codes may be constructed with the help 
of beautiful properties of GHZ states. 

There still remain many interesting problems about our work in this paper. 1) 
Devise a theoretical framework for describing quantum error correction codes in the 
Hamiltonian formulation without involving the stabilizer formalism. Our typical ex- 
amples are mainly based on algebraic properties of GHZ states. 2) Study physical 
properties of our Hamiltonian models and explore their applications to quantum error 
correction. As a comparison, quantum error correction code in the Hamiltonian model 
based on the stabilizer formalism is its degenerate ground state, and there is a gap be- 
tween this ground state and its first excited state. 3) Study experimental realizations of 
these Hamiltonian models on optical lattices or atomic spin lattices [681 [69] . The uni- 
tary basis transformation from the product basis to GHZ states can be experimentally 
performed in view of the work [JI] . 4) Study applications of these Hamiltonian models 
to develop a theoretical framework for integrable quantum computation or topological- 
like quantum computation. Note that the unitary evolution operator in the Hamiltonian 

2 The author thanks Zohar Nussinov for the email correspondence on this point. 
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model for a quantum error correction code using GHZ states can be chosen as a solution 
of the quantum Yang-Baxter equation. 5) About quantum error correction codes using 
GHZ states, there exist at least three types of questions to be answered: devise a gen- 
eral framework for yielding nonadditive codes; study quantum circuits for encoding, 
error correction and decoding; explore fault-tolerant quantum computation with the 
help of properties of GHZ states, for example, channel-adapted fault-tolerant quantum 
computation for the amplitude damping channel. 
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